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Abstract  -  The  existing  methods  developed  in  URTURIP 
Technique  (Ultrasound  Reflection-mode  Tomography  Using 
Radial  Image  Processing)  can  be  considered  as  a  part  of  a 
reconstruction  family  called  Maximized  Weighted  Averages 
Technique,  and  it  is  the  object  of  this  paper.  It  is  shown  that  a 
reconstruction  can  be  viewed  as  a  filter  and  that  the  three 
commonly  used  methods  can  be  rewritten  with  the  same 
formulation.  Moreover,  interesting  new  methods  can  be 
developed  by  using  other  adapted  filters. 
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I.  Introduction 

As  a  part  of  our  project  concerned  with  the  development 
of  an  ultrasound  scanner  dedicated  to  limb  study,  the 
URTURIP  Technique  (Ultrasound  Reflection-mode 
Tomography  Using  Radial  Image  Processing)  has  been 
developed  [1].  It  consists  in  using  classical  B-scan  images 
instead  of  projections  [2-5]  and  gives  qualitative  images 
instead  of  quantitative  images.  The  final  goal  of  this  project  is 
the  2D  and  3D  reconstruction  of  anatomical  structures  at  limb 
level  by  using  echographic  image  processing.  The  developed 
process  consists  of  several  successive  steps  like  :  multiple 
reflection  removing  [6],  2D  reconstruction  [1],  segmentation 
[7],  contour  association  [8],  contour  interpolation  [9],  3D 
reconstruction  and  visualization  [10].  It  has  been  validated  by 
in  vitro  experiments  on  anatomical  pieces  of  limbs  of  new¬ 
borns  using  a  simple  acquisition  system  prototype  [10].  Now, 
in  vivo  experiments  are  carried  out  thanks  to  a  calibration 
technique  [11]  which  allows  a  good  positioning  of  the  probe 
and  an  accurate  determination  of  the  rotation  centre. 

Most  of  the  existing  methods  developed  in  URTURIP 
Technique  can  be  considered  as  a  part  of  a  reconstruction 
family  called  Maximized  Weighted  Averages  Technique,  and 
it  is  the  object  of  this  paper.  It  is  shown  that  a  reconstruction 
can  be  viewed  as  a  filter  and  that  the  three  commonly  used 
methods  can  be  rewritten  with  the  same  formulation. 
Moreover,  several  new  interesting  methods  can  be  developed 
by  using  other  filters.  It  allows  adapting  the  reconstruction  to 
an  application  and  choosing  a  good  contrast  between  the 
different  structures  to  be  observed. 


II.  The  urturip  technique 

A.  Principle  of  reconstruction 

The  Ultrasound  Reflection-mode  Tomography  Using 
Radial  Image  Processing  (URTURIP  Technique)  principle  is 
to  utilize  radial  B-scan  images  instead  of  projections  as  most 
other  methods  do  [2-5].  In  comparison,  fewer  radial 
directions  are  needed,  it  is  less  time-consuming,  but 
qualitative  images  instead  of  quantitative  ones  are  computed. 

Let  i=l,..,N  be  N  radial  images  obtained  from  N 
angulary  equidistant  directions  around  the  rotation  center 
where  Lj  (k,  1)  denotes  the  luminance  of  the  (k,  l)-pixel  on 
the  image. 

An  adjustment  step  consists  in  turning  each  Lj  image 
around  the  rotation  center  with  its  own  acquisition  angle,  to 
construct  N  adjusted  images  in  such  a  way  that  a  pixel  (x, 
y)  on  each  image  corresponds  to  the  same  real  point  of  the 
cross-section. 

Then,  a  method  based  on  the  URTURIP  Technique 
reconstructs  an  image  L  by  a  combination  of  the  N  adjusted 
images  ,  i.e.  L(x,  y)  =  f(Li(x,y))  ,  i=l..N,  where  f  denotes 
the  reconstruction  method. 

B.  The  commonly  used  methods 

The  three  commonly  used  methods  are  the  method  of 
maxima,  the  method  of  averages  and  the  method  of 
maximized  averages  [1].  They  are  defined  as  follows. 

The  method  of  maxima  is  given  by  : 

V(x,y)  L(x,y)  =  maxL^xj)  i  =  l..N 

i 

The  method  of  averages  is  given  by  : 

i  N 

V(x,  y)  L(x,y)  =— y^L;(x,y) 

i=l 

The  maximized  averages  method  combines  the  respective 
advantages  of  previous  two  basic  method  in  order  to  obtain  a 
good  image  quality  even  with  few  radial  images.  Its  formula 
is  as  follows: 

V(x,y)  L(x,y)  =  maxLj(x,y)  i  =  l..N,  N>5 
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with 


1  ni+1 

L'i(x’y)=^— [ELk©(x’y)  i=1- 


N 


n,  =  l  -  m 


ni+1  -  i  +  m 


m  =  < 


J=ni 


N  , 

- 1  otherwise 

4 


k(j)  = 


j  if  je[l..N] 
j-N  if  j  >  N 
j  +  N  if  j<l 


where  E(z)  denotes  the  integer  part  of  z. 

III.  The  maximized  weighted  averages 
A.  Formulation 

The  three  previous  method  can  be  described  by  a  same 
formulation  defined  as  follows  : 


V  (x,y)  e  L  L(x,y)  =  Max  Lj(x,y) 


1  =  1.. N 


1  N 

with  L'i(x,y)  =  — ^  C(n(j))Lj(x,y) 


where 


"i  j) = 


Ec(j)  i'1 

j=-N/2+l 


j-i  +  N  if  j<Sup|i-y,  1 


J-i 


if  Sup|i-y,  lj<j<lnf[i+y,  N 


N 


j-i-N  if  j  >  Inf|  i  N 


where 


n(j)  = 


IH 

if 

j-i 

1  <  — 

1  2 

N  |  j  i| 

if 

|j-i| 

N 

>  — 
2 

B.  Rewritten  of  the  existing  methods 

It  clearly  appears  that  the  three  commonly  used  methods 
correspond  to  three  common  filters  (uniform,  impulsive  and 
square  window).  They  can  be  rewritten  with  the  same 
formulation  and  the  corresponding  filter  C  or  C  which 
defines  the  method  used. 

The  method  of  averages  (fig.  1)  is  defined  by  : 

C(k)  =  1  V  k  =  -N/2  +  l,..,N/2  or 

C*(k)  =  1  V  k  =  0, . . , N  / 2 

The  method  of  maxima  (fig.  2)  is  defined  by  : 
fl  if  k  =  0 

C(k) = \  r  ,  r  ,  or 

'0  ke[-N/2  +  l,N/2]-{0} 

1  if  k  =  0 
0  if  1  <  k  <  N/2 


C  (k)  = 


The  method  of  maximized  averages  (fig.  3)  is  defined  by  : 
fl  if  -  m  <  k  <  +m 

C(k)  =  \  r  r  i  or 

7  [0  ke[-N/2  +  l,N/2]-{0} 


C  (k)  = 


1  if  -  m  <  k  <  +m 
0  if  1  <  k  <  N/2 


-N/2 


N/2 


C  is  a  filter  defining  each  method.  Moreover,  C  being  a 
symmetric  filter,  its  restriction  C  on  [0,  N/2]  can  be 
considered  and  the  formulation  becomes  : 


V  (x,y)  e  L  L(x,  y)  =  Max  Li(x,y) 


i  =  1..N 


N 


with  Lj(x,y)  =  - 


N/2 


^C*(n(j))Lj(x,y) 


c*(o)+2y;c*(j)  j=i 


N/2  -1  N/2 


Fig.  1  :  C  and  C*  for  the  method  of  averages 
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Fig.  2  :  C  and  C*  for  the  method  of  maxima 


images.  It  can  be  adapted  to  different  applications  and  depth 
field  configurations. 
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Fig.  3  :  C  and  C*  for  the  method  of  maximized  averages 


Fig.  4  :  C  and  C*  for  a  exponential  filter 
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C.  New  methods 

Of  course,  in  this  family  of  Maximized  Weighted  Averages 
Methods ,  other  methods  can  be  developed.  In  the  context  of  a 
comparative  study,  different  methods  have  been  tested  using 
different  filters  [12]  like  Haming,  Hanning,  median,  ... 

The  Maximized  Weighted  Averages  Methods  using  an 
exponential  filter  (fig.  4)  or  a  gaussian  filter  (fig.  5)  give  quite 
interesting  and  good  results  and  are  currently  used  in  our 
project.  They  are  defined  by  : 

C(k)  =  e_0lkl  for  k  e  [-  N  /  2  +  1,N  /  2]  or 

C*(k)  =  e_ok  forke[0,N/2]  and 

C(k)  =  e“k2/2<j2  for  ke[-N/2  +  l,N/2]  or 

C*(k)  =  e-k/2°  for  k  e  [0,N  /  2] ,  where  g>0  is  an 
interesting  coefficient  to  weight  the  influence  of  radial 


Fig.  5  :  C  and  C*  for  a  gaussian  filter 


(a)  gaussian  filter  (b)  exponential  filter 

Fig.  6  :  2D  reconstruction  at  thigh  level 


Fig.  6  shows  a  reconstructed  cross-section  at  thigh  level 
with  maximized  weighted  averages  methods  respectively 
based  on  a  gaussian  filter  and  an  exponential  filter.  It  gives 
good  contrast  to  distinguish  the  different  anatomical 
structures. 
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IV.  Conclusion  [12]  L.  Chazallon,  Calibration  du  systeme  d’acquisition  d’un 

scanner  a  ultrasons  medical  base  sur  la  technique  URTURIP, 
In  the  context  of  a  project  of  the  development  of  an  Rapport  de  DEA  SIBM,  Universite  d’ Angers,  (1998). 
ultrasound  scanner  dedicated  to  limb  study,  the  URTURIP 
Technique  has  been  developed. 

In  this  paper,  Maximized  Weighted  Averages  methods  are 
described.  Thanks  to  a  unique  formulation  all  the  three 
commonly  used  methods  can  be  rewritten.  It  is  shown  that 
they  correspond  to  different  classical  filters.  Moreover,  in  the 
same  category  of  Maximized  Weighted  Averages,  new 
interesting  methods  can  be  developed  by  using  other  filters. 
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